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Abstract 

A first-principles-derived approach is developed to study the effects of un- 
compensated depolarizing electric fields on the properties of Pb(Zr,Ti)03 ul- 
trathin films for different mechanical boundary conditions. A rich variety of 
ferroelectric phases and polarization patterns is found, depending on the in- 
terplay between strain and amount of screening of surface charges. Examples 
include triclinic phases, monoclinic states with in-plane and/or out-of-plane 
components of the polarization, homogeneous and inhomogeneous tetragonal 
states, as well as, peculiar laminar nanodomains. 
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Ferroelectric thin films are of increasing technological interest because of the need in 
miniaturization of devices [1]. An intriguing problem in these films concerns their polar- 
ization patterns. For instance, various patterns have been recently proposed from mea- 
surements, phenomenological theory and first principles. Examples are out-of-plane mon- 
odomains [2-5], 180° out-of-plane stripe domains [5,6], 180° and 90° multidomains that are 
oriented parallel to the film [7], and microscopieally-paraelectric phases [4]. The fact that 
dramatically different patterns have been reported for similar mechanical boundary con- 
ditions supports a concept discussed in Refs [4,8,9], namely that they arise from different 
electrical boundary conditions. More precisely, real thin films are likely neither in ideal 
open-circuit conditions — for which unscreened polarization-induced surface charges can 
generate a large depolarizing electric field along the growth direction [10,11] — nor in ideal 
short-circuit conditions — that are associated with a vanishing internal field resulting from 
the full screening of surface charges — , but rather experience a situation in between. As a 
matter of fact, a reactive atmosphere can lead to a partial compensation of surface charges 
in films with nominal ideal open-circuit conditions [5] , while metallic or semiconductor elec- 
trodes "sandwiching" films do not always provide ideal short-circuit conditions — resulting 
in a (non-zero) internal field, whose magnitude depends on the element from which these 
electrodes are made [4]. The amount of surface charges' screening in thin films can thus vary 
from one experimental set-up to another, possibly generating different polarization patterns. 

A precise correlation between this amount of screening and the morphology of the polar- 
ization pattern, and how this correlation depends on mechanical boundary conditions, are 
still lacking nowadays despite their obvious importance. Similarly, atomic-scale details of 
multidomains - and their formation mechanism - are mostly unknown in ferroelectric thin 
films. One may also wonder if some uncompensated depolarizing fields can yield unusual 
polarization pattern, such as monodomains made of ferroelectric phases that do not exist in 
the corresponding bulk material. Candidates for these latter anomalous features to occur 
are films made of perovskite alloys having a composition lying near their morphotropic phase 
boundary (MPB), because of the easiness of rotating their spontaneous polarization [12-15]. 

In this Letter, we develop a first-principles-based scheme to investigate the effects of 
uncompensated depolarizing fields on the properties of Pb(Zr,Ti)03 films near their MPB, 
for different mechanical boundary conditions. Answers to the problems summarized above 
are provided. In particular, we find a rich variety of ferroelectric phases, including unusual 
triclinic and monoclinic states. We also observe complex nanodomains, and reveal their 
formation and atomic characteristics. 

Specifically, our Pb(Zri_j.Tia.)03 (PZT) thin films (i) are grown along the [001] direction 
(to be chosen along the z-axis); (ii) are "sandwiched" between non-polar systems (mimicking, 
e.g., air, vacuum, electrodes and/or non- ferroelectric substrates); (iii) have Pb-0 terminated 
surfaces; and (iv) have a 50% overall Ti composition. Such low-dimensional structures are 
modeled by large periodic supercells that are elongated along the z-direction, and that 
contain a few number of B-layers to be denoted by m — with the alloyed atoms being 
randomly distributed inside each of these m planes. (The resulting films thus have a thickness 
^ 4m A). Typically, we use 10 x 10 x 40 periodic supercells with m around 5. The non-polar 
regions outside the film are thus altogether 40 — m lattice constant thick along the growth 
direction, which allows well-converged results for the films properties. The total energy of 
such supercell is used in Monte-Carlo simulations, and is written as: 
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Etot{{Ui}, {Vj}, 77, {ai}) = EHeffii^i}, {Vj}, 77, {(7i}) 

- Z!/^27r-^ < Mj> >3 Mi,^ (1) 
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where Eneff is the (aUoy effective Hamiltonian) energy of the ferroelectric film per se. Its 
expression and first-principles-derived parameters are those given in Ref. [14] for bulk PZT. 
Uj are the local soft modes in unit cells i of the PZT film — which are directly proportional 
to the electrical polarization and whose components along the z-axis are denoted as Ui^z- 
{vi} are inhomogeneous strain-related variables inside these films [16], while {r}} is the 
homogeneous strain tensor. As indicated in Ref. [17], the form of {rj} (in Voigt notation) is 
relevant to two cases of interest, namely stress-free vs. epitaxially strained (001) films. In the 
former case, all the components of {77} fully relax. On the other hand, the second situation is 
associated with the freezing of three components of {rj} — i.e., tiq — and rji — r)2 = S, with 
5 being the strain resulting from the lattice mismatch between the film and the substrate 
- while the other components relax during the simulations [17]. {cTj} characterizes the 
atomic configuration, i.e. (Tj=+1 or —1 corresponds to the presence of a Zr or Ti atom, 
respectively, at lattice site i in the PZT film. The local modes and the inhomogeneous 
strain-related variables are forced to vanish in the supercell areas located outside the PZT 
films. Therefore, a depolarizing field implicitly occurs inside the film if this latter has a 
component of its polarization along the growth direction. The second term of Eq (1) mimics 
the effects of an internal electric field — that arises from the partial or full screening of 
polarization-induced charges at the surfaces, and that is thus opposite in direction to the 
unscreened depolarizing field — on the films properties. This energetic term is dependent 
on the Z Born effective charge [16], the a lattice constant, the e^o optical dielectric constant 
of Pb(Zro.5Tio.5)03) ^-nd the average of the z-component of the local modes centered at the 
surfaces (denoted by < Uj^z >s, and that is self-consistently updated during the simulations) 
[18]. This second term is also directly proportional to a /3 parameter that characterizes 
the strength of the total electric field inside the film. Specifically, P — corresponds 
to ideal open-circuit boundary conditions with E,^ having its maximum magnitude (when 
polarizations lie along the z-axis), while an increase in f3 lowers this magnitude. 

The value of f3 resulting in a vanishing total internal electric field is dependent on the 
supercell geometry, and in particular on the number of its non-polar layers [10]. This ideal 
short-circuit (3 is denoted as Psc in the following. It is found to be 0.69, for stress-free films 
associated with m — 7 and a supercell periodicity of 40 lattice constant along the z-axis, 
by fitting the T = 10 K predictions delivered by Eq (1) to a single result of Ref. [3] (that 
constructed a microscopic Effective Hamiltonian for thin films under stress-free and ideal 
short-circuit boundary conditions); that is, the polarization in the film layer that is further 
away from the surfaces is along the z-axis and has a magnitude equal to the one in the bulk 
[19]. (Note that Pb(Zro.5Tio.5)03 bulk is tetragonal with an average local mode of 0.1072, in 
lattice constant units, at T=10 K [14]). We found that our resulting T = 10 layer-by-layer 
profiles of the local modes for stress-free ultrathin films, that are under ideal short-circuit 
electrical boundary conditions, are remarkably similar to those of Ref [3]. For instance, the 
polarization at the surfaces is significantly enhanced with respect to the bulk, and increases 
as the film thickness decreases. Ref. [3] and our proposed scheme also both predict that the 
layers exhibiting the smallest polarization are those next to the surface layers, and that the 
layer that is further away from the surfaces has a larger polarization for ttt, = 5 than for 
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m = 3 or 7 (thus exhibiting a non-monotonic behavior versus thickness). 

Having demonstrating that our approach can capture subtle details, we now apply it to 
study the effects of uncompensated depoienizing fields on the ground-state of Pb(Zro.5Tio.5)03 
ultrathin films. The thickness is kept fixed at m = 5, while different mechanical boundary 
conditions are adopted: (a) stress-free (Fig la); (b) a tensile strain S — -1-2.65%, cor- 
responding to an increase of 0.02 Bohr for the in-plane lattice constant with respect to 
Pb(Zro.5Tio.5)03 bulk (Fig lb); (c) a compressive strain d = —2.65% (Fig Ic). More pre- 
cisely. Figs 1 show the predicted (absolute value of the) Cartesian components < >, 

< Uy > and < Uz > — along the [100], [010] and [001] directions, respectively — of the 
average of the local mode vectors in the film, as a function of P/Psc and at T=10K. Figs 1 
also display the behavior of um, which is defined as um — ■sj< ul + Uy + > and thus 
provides a measure of the local polarizations. 

Under stress-free conditions, the film has a spontaneous polarization aligned along the 
z-axis for (large) values of (3 that correspond to a screening of at least 98% of the polarization- 
induced surface charges. This results in a tetragonal state to be denoted by T^.. On the other 
hand, when P becomes smaller than ~ O.dOAPsc, the internal field along the growth direction 
would be too strong to allow an out-of-plane component of the local mode. As a result, the 
polarization aligns along an in-plane < 010 > direction. The corresponding ferroelectric 
phase is denoted as T^^. The most striking result for stress-free PZT films is the polarization 
path when going from T^. to Ty. As P/Psc decreases from 98% to 90.4%, the polarization 
continuously rotates and passes through three low-symmetry phases: a so-called monoclinic 
state [20] — occurring for 0.932 < P/f3sc < 0.98, and for which < Uy > and < Ux > are 
nonzero, equal to each other and smaller than < Uz >', & triclinic Tr phase, for {3/ f3sc ranging 
between 92.2 and 93.2%, that is characterized by a local mode with non-zero and different 
Cartesian components; and a so-called monoclinic ground-state [20] , when [3 / Psc ranges 
between 90.4 and 92.2%, for which < > vanishes while < Uy > becomes larger than 

< Uz >■ Interestingly, neither the Tr nor the phase exists in the temperature- versus- 
composition phase diagram of PZT bulk! Decreasing dimensionality thus leads to unusual 
phases because of residual depolarizing fields. Moreover, all the low-symmetry states of 
Fig la are only stable for compositions lying near the MPB of bulk PZT since we found that 
the M^, Tr and M^ phases disappear in favor of T^ for higher Ti concentration. 

Comparing Fig. lb with Fig. la reveals that films under a tensile strain react to depo- 
larizing fields in a dramatically different way than stress-free films. In particular, the end- 
member phases (T^ and T^) of Fig la both disappear. The reason for the vanishing of the 
(ideal-short-circuit-derived) T^ phase is that tensile strains favor in-plane components of the 
local modes, because of the well-known coupling between strain and polarization [16,17,21]. 
Consequently, T^ transforms into a M^ phase for large enough r]i = ri2 = S strains. These 
epitaxial constraints — and more precisely the fact that 772 can not be different than r/i — 
generate an in-plane component of the polarization along the [100] direction in addition to 
a larger component along the [010] direction, for small j3. The (ideal-open-circuit-derived) 
Tj, state thus becomes a Mc phase when going from stress-free to tensile conditions [22] . 

Conversely, a large enough compressive strain annihilates the (in-plane) < > and 

< Uy > components of the local mode for any (3 (see Fig. Ic). Two macroscopically-diiiereni 
phases result from this annihilation: a paraelectric P phase for /? smaller than 0.822(3sc 
and a ferroelectric tetragonal T^ phase for larger (3. T^ can be further separated into two 
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microscopically-diSerent phases. For (5/(5sc > 0.884, the local polarizations all point along 
the growth direction and have similar magnitude, since < > is nearly equal to Um- The 
resulting state is referred to as Tj''^. On the other hand, when (5 ranges between 82.2% and 
88.4% of I3sc-, < > becomes smaller than um- This characterizes a locally- inhomogeneous 
polar state to be denoted by TJ'^. The transition between T'j'*) and TJ'^ is of first-order, as 
revealed by the sudden jump of both < Uz > and um- Another striking feature that can 
be extracted from Fig. Ic when looking at < ^2 > and um is that the macroscopically- 
paraelectric P phase has a relatively large magnitude for its local polarizations! 

In fact. Figs 1 tell us that the only phases for which the magnitude of the macroscopic 
polarization significantly differs from that of the local polarizations are TJ'^ and P. Un- 
der thermal equilibrium, complex polarization patterns (if any) should thus only exist for 
compressively-strained ultrathin films, while large ferroelectric monodomains are likely to oc- 
cur in ultrathin films that are under stress-free or tensile mechanical boundary conditions. 
In order to provide a detailed microscopic insight of compressively-strained thin films. Figs 2 
(a), (b) and (c) display a snapshot of (very large) 24 x 24 x 40 supercell simulations yielding 
a T^^\ Tj*^ and P phase, respectively. Fig 2(a) confirms that T^^^ is locally homogenous. 
Tj^) is likely the phase observed in Ref [2] since this latter also occurs for compressive 
strain and large screening of the surface charges. More striking. Fig 2(b) reveals that Tj*) is 
characterized by the formation of nanodomains having local dipoles that are aligned in an 
opposite direction with respect to the macroscopic polarization. These nanodomains prop- 
agate throughout the entire thickness of the film but are laterally confined. We believe that 
Tj*) is the phase experimentally seen in Ref [5], and denoted there as F^, since this latter 
exhibits features that are consistent with our inhomogeneous Tj*\ namely: (i) compressive 
strain conditions and partially compensated depolarizing fields; (ii) a non-zero spontaneous 
polarization along the z-axis; and (iii) broad diffraction peaks. In the P state (see Fig 2c), 
the nanodomains have laterally "percolated" , resulting in the formation of nanoscale 180" 
(out-of-plane) stripe domains. Such peculiar multidomains are consistent with the so-called 
Fq, and F/3 phases observed in Ref [5]. Another interesting feature revealed by Fig. 2(c) 
— which is, to the best of our knowledge, the first atomic-scale picture of nanodomains in 
ferroelectric thin films — is that each domain is terminated at one surface with local dipoles 
aligned in opposite in-plane directions, while the other surface has out-of-plane dipoles. This 
contrasts with a commonly accepted picture of out-of-plane 180° domains in magnetic films 
for which both surfaces have in-plane magnetizations [23]. Finally, Fig. 2c also shows that 
the predicted period A of the laminar domains is ^ 8aj„, where a^^ is the in-plane lattice 
constant. This remarkably agrees with the measurements of Ref. [5] yielding A = ?,1A for 
compressively-strained 20A-thick PbTiOa films. (Other calculations we performed — using 
n X n X 40 supercells with n=10, 16, 18, 20 and 24 — confirm that A ~ 8aj„ or 9ai„) 

We thank D.D. Fong, R. Haumont, J. Junquera, 1. Naumov, S.K. Streiffer and C. Thom- 
son, for having stimulating our interest in the present topic. This work is supported by ONR 
grants N 00014-01-1-0365 and 00014-01-1-0600 and NSF grant DMR-9983678. 
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FIGURES 



FIG. 1. Cartesian components of the film average of the local- mode vectors, as a function of the 

(3/Psc parameter in (001) Pb(Zro. 5X19.5)03 ultrathin films having a m = 5 thickness, at T = lOK. 
Part (a) displays the predictions for stress-film films. Parts (b) and (c) show the results for a 
tensile and compressive strain of 2.65%, respectively. The arrows indicate the values of /5//?5c 
around which phase transitions occur. The results for (i/ f3sc down to zero are not shown since 
they are identical to those associated with the lowest values of P/Psc displayed here. 10 x 10 x 40 
super cells are used. 

FIG. 2. Three-dimensional polarization patterns in (001) Pb(Zro.5Tio.5)03 films having am = 5 
thickness and that are under a compressive strain of —2.65%, at T = lOK. Part (a), (b) and (c) 
corresponds to a p/Psc parameter of 94.5% {T^'^^ phase), 87.7% (tJ*^ phase) and 80.8% (P phase), 
respectively. The bottom of Fig. 2c shows the projection of the 3D picture into a X-Z plane. 
The red (respectively, blue) color indicates local dipoles having a positive (respectively, negative) 
component along the z-axis. 24 x 24 x 40 supercells are used. 
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This figure "Fig2.png" is available in "png" format from: 
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